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Abstract 
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calculating three undetermined functions which defined a solution in the 
LTB model are presented. One example of exact nonhomogeneous model is 
studied. The limit transformation to the FRW model is shown. 
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1 The Introduction 

The LTB model is one of the most known spherical symmetry model in general relativity. It was 
created by Lemaitre jy, Tolman and Bondy || during the period of time from 1933 to 1947. 
The exact solution have been obtained by Bonnor Q| in 1972 and Q in 1974. The LTB model 
represented one of the simplest nonhomogeneous nonstationary cosmological models and due to this 
fact is used to study some new ideas in the cosmology. 

The present interest to LTB was risen by the observations shown the fractal structure of the 
Universe in the large scale @, 0- The review of the topics in Cosmology is presented in In 
a set of papers the LTB model is used to study the observational datas and redshift as a main 
cosmological test || - O ]. 

The central problem of using the LTB model is in calculation of three undetermined functions 
which defined the solution. There is a set of ways how to solve this problem from the physical point 
of view in the mentioned papers. This article is devoted to the mathematical point of view on this 
matter. 

2 The LTB Model 

This section is devoted to presentation the Lemaitre- Tolman-Bondy model and section 2 of the 
paper Q is cited. The co-moving system of coordinate is used in this model where the interval has 
the form 

ds 2 {r, t) = -e x ^dr 2 - e w(r ' 4) (d9 2 + sin 2 ed<£) + dt 2 . (1) 

A(r, t) and u>(r, t) are metrical functions j|] definding the solution. In the co-moving system of 
coordinate with the line element (Q) the energy-momentum tensor 

„ R dx a dx? 

has only one non zero component 

T 4 4 = p Tp — 0, a or (3 = 4. 
Using them together with Dingle results [ll| we obtain the system of equations of the LTB model: 

%- K T} = e- UJ -e- x — +d) + -c0 2 -A = (3) 
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8nTi = e-" - e~ x U + h 2 ^ ) + — + '— - A = Ntt/; Ci 
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where 
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The equation (o) has the solution 



(7) 

e X = e^-^, (8) 
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where f(r) is undetermined function. Substituting (^j) into (||) we obtain 



e"(o) + ^ 2 -A) + [l-/ 2 (r)] =0. 



This equation is integrated twice. First integral gives the equation 

,2 



e 3 " /2 (f - - |a) + 2e^ 2 [1 - / 2 (r)] = F(r), (10) 



and the second one gives the equation 

de"/ 2 



^/f 2 (r)-l + ii^(r)e— / 2 + fe^ 



= t + F(r) (11) 



The equations ( |To| ) and ( |Tl] ) hold undetermined functions F(r) F(r). The substitution of (|J) into 
(0) together with ( |i~0| ) gives the equation for density 

3 The Cauchy Problem for the LTB Model 

Before we study the LTB model let us introduce the follow characteristic values: a velocity of light 
c, an observational meaning of the Habble constant H, a characteristic time 1/H and characteristic 
length c/H. We use the co-moving system of coordinates in the LTB model, so the radial coordinate 
r has a sense of Lagrangian mass coordinate pq,p7y. Two dimensionless variables \l and r are 
dcfinded by the rules 

fi=— T = Ht, (13) 

m 

where m is a full mass of the "gas". The dimensionless Habble function h(jJ,, t) and density 5(/x,t) 
will be also used: 

i / \ H(p,t) xt s p(r,t) 

where p = p(0, 0), H(0, 0) = ff. 
Let us write the interval (|l|) as 

ds 2 (r, t) = -Ae x{r > t] dr 2 - Be u(r ' t) (d6 2 + sin 2 0d<f>) + c 2 dt 2 , (15) 

where two constants A and B are introduced to take into account the fact that ( |l5| ) is dimension 
equation. 

The dimension of [ds 2 ~\ is L 2 , dimension of [A] is L 2 M~ 2 and dimension of [B] is L 2 , so 

The interval (||) has now the form 

C ^ ds 2 (r, t) = -e^W 2 - (dfl 2 + Sl n 2 9d^ + dr 2 , (17) 



Together with metrical functions w(/Lt, r) and A(/x,r), introduced by Tolman, it is conveniently to 
use the Bonnor's function Q 

i?( M , T ) = e "(/vr)/2_ (18) 
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In the Bonnor's notation the interval (17) takes the form 

^ ds\ix, t)=- [R '^ 2 T)]2 dfi 2 - R 2 (fi, r) (d9 2 + sin 2 6dcj>) + dr 2 (19) 

As it is shown in |ll| and fj"7| , the Bonnor's coordinate R{p., r) has a sense of Euler coordinate, so 
the equation ([lj) correlates geometrical radius of the sphere R([i, r) where the particle is located, 
and the Lagrangian coordinate fi of this sphere. 

To describe the radial motion we will use the Habble function connected with variation of the 
length dl: 

h =i < 20 » 

where, according to the (p"7|), for dl 2 we read: 

dl 2 = e x ^dn 2 =(jdf?j . (21) 

By the substitution ( |2l| ) into the definition of the Habble function (|2^) we obtain 

, , , A(m,t) R! d\nR' 

M^) = ^ = - = — (22) 



By the integration of the equation (22) we obtain the formula for metrical function A(/i,r): 

T 

A(m, t) = 2 / h{p, T)dr + A(m, 0) (23) 



A solution in the LTB model is defined by the functions f(r), F(r) and P(r). These functions 
are obtained in the process of solution of the system of PDE, so to define them the initial/boundary 
conditions should definitely be used. The metrical function cj(/i,t) is the solution of the equation 
( pT| ) . The equations of the LTB model are obtained in Q and solved in the parametric form for the 
three cases f 2 (n) < 1, flitt) = 1 and / 2 (m) > 1 m and [||- 

The equations @ - ([ll] ) are valid for every r, and due to this fact in the Cauchy problem they 
definde the functions / 2 (/i), F(fi) and F(/i) at the moment of time r = 0: 



f 2 ( M ) - 1 = e u »W + ^ 2 (fx) a) 



(24) 



F(/i) = c 3 ^/ 2 - Ha) + 2e^>/ 2 [1 - f 2 ( M )] (25) 

F(M)= / I (26) 
c . /o, /2 V/ 2 (M) - 1 + lF(ji)e-^)/2 + AgOoM 

At the time r = the equation (^) defindes the function Ao(/z): 
Substituting ( p4| ) into (p5|), we obtain 



(28) 
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Comparing (§) - (0) with (||) - ©, we find out that 

f( f ,) 1 = e-M ^o(M) + ^o 2 (i») - a) 
e^ r) foi(A*,r) + |a; 2 (/i,r)-A) 



(29) 



and 



(30) 

are not dependent on time. Let's use the previous results to calculate the functions F(/i) and integral 
in the equation (fLy). Substituting the definitions (p3) and (pq) into (pq) we obtain: 



F(/,) = ± / (31) 

w (0) 

The function F(/i) is equal to zero at the moment of time r = according the definition. 
Substituting the right part of the equations ( p^ ) and (^) into the (^6|) , we obtain the equation 

± / — =± / — +r. (32) 

J CO J 10 

U)(ft,0) w (0) 

This analysis of the LTB model shows that the functions 

w (^T)lr=o = <*>o(m) w(/i,r)| T=0 = c0 o (^) 

&O,t)| t=0 =w (m) ! 
and constants 

w(m,0)| m=o =w o (0) w( M ,0)| M=o =d;o(0) ) 

(34) 



(33) 



w(M, 0)| p=0 =u>o(0), A, 

are included into the definitions ( pi[ ) - ( p6| ) and they form the initial conditions of the Cauchy 
problem for the equations (||) - (|j). In accordance with (p7|) the function Ao(/i) is not include in the 
set of initial conditions. 

Substituting into the ( [l2"| ) , we obtain the general expression for the density of " gas" in the 
LTB model: 



3[u)o{n)]' 



•• r s 1.9/ \ A 



87rp(^, r) = r x 

- 2 [woO*)]' - 2wo0i) [wo(m)]'| (35) 



The function o;(/x, r) from the equation (|3^) is the solution of the equation gj. 

4 The / 2 = 1, A = solution 

The solution of the main equation of the LTB model (|ll]) will be obtained in this section for one 
special case 

f = 1, A = 0. (36) 
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This case allows simple analytical solution. The law of the density, Habble function, and cosmological 
parameter will be also studied here. 

The condition / 2 (/i) = 1, as it goes from the (|34|), correlates the initial conditions as follows: 

wo fa) + - 0. (37) 

Due to this reason, the general number of the initial condition (|3~i| ) is decreased by one unit. To 
obtain the solution we'll start with the solution of the equation (|10|). ^From it it's obvious that: 



e<". (38) 



In the general case the function F is defined by the equation (p5|). Whit the help of ([56j) it takes 
the form 

F = ^e 3 "°/ 2 . (39) 
The equation (^) with conditions (|3^) takes the form 

a„3u;/4 Q 

C -o/ 4 ^_ = ± ^ . (40) 

Integrating (|4(]) we obtain 

e f[™]=±^ or + F(/x). (41) 
We find out the function F from the initial conditions then time r = 0: 

F = 1. (42) 

The solution now is: 

e j[w-wo] = i±^ QT . (43) 
In the Bonnor's notation the equation (0) becomes as follows: 

where the upper sign correspond to the expansion into the infinity from the initial condition and 
lower sign correspond to the collapse from one. The solution with "-" describes time of collapse 
depening on the initial mass coordinate: if the particle has a Euler coordinate Ro(fJ-) at the moment 
t = 0, the time of collapse is equal to 

f( M ,0) = — ^- (45) 
We obtain the density corresponding to this solution by the substitution ( fi"3| ) into fl35|). Let's denote 

(46) 

87r5(/i,T) = - -2 |t= , (47) 



(1 ± f w r) 1 ± 



1 ± |w r 

Using the definition (|2^) we find out the Habble's function for the solution ([i"3|): 

T) = ± ^° + \ =p i • -^L_ (48) 



(> 



Knowing the law of the density and Habble's function we obtain now the formulum for the cosmo- 
logical parameter. By the definition 

0=-, where p c =^, (49) 
p c oitG 

where H and p mean the Habble function and the density at the moment of the observation. Let's 
assume that this moment is r = 0. The density, the critical density and Habble function at this 
moment are: 

8tt5(m,0) =w (j^o + ^j , (50) 
Mft0) = ^ 2 ( ft 0). (51) 
%,0) = ±Qw + ^ (52) 
Q (/i ,0) = gg% fo> + " (53 ) 
These functions dependent on time as follow: 



* c(/i ' T) = ^§p~ h2 ^ T) - (54) 



16Gp . /3 . 

II = = OJq -UJn + V 

3H 2 V4 



(1 ± ftior) [1 ± (|«o + f) r] 



(55) 
(56) 



{±4 (1 ± |c0 o r) (|w + T [l ± (|w + i/) r] } 

5 FRW model 

Every nonhomogeneous solution of the LTB model must include the FRW model as a limited 
case when the density and Habble function are not dependent on space coordinate for all moments 
of time. Let's study in which case the present solution is reduced to the FRW model? Only one 
condition satisfied this request: 

ujq = const so, vq(p) = 0. (57) 



It goes from mw and (57) that in this case 



^S=U-^) , (58) 



4 V 1 ± jujot 

^ = 4t-¥- (59) 

2 1 ± fuj Q T 

^From these equations it goes: 

8ttS = 3/i 2 , (60) 

We obtain the meaning of the loq from (|58| ) : 

u; Q = ±2^. (61) 

v = together with ( |6l| ) give 

n = ^. (62) 
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